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The resonant generation of relativistic plasma waves and plasma wave guiding by two co- 
propagating laser pulses has been studied. By proper timing between guiding and driver pulses, 
a resonant interaction occurs, which generates a high-amplitude plasma wave over a longer length 
than is possible with either of the laser pulses individually. The growth of the plasma wave is inferred 
by the measurement of monoenergetic electron beams with low divergence that are not measured 
by using either of the pulses individually. This scheme can be easily implemented, and allows more 
control of the interaction than is available to the single pulse scheme. 

The laser wakefield accelerator (LWFA) |IJ consists of a high-intensity laser pulse that generates an electron plasma 
wave with relativistic phase velocity, and oscillating at the electron plasma frequency , u p = y / e 2 n e /m e eo, where n e 
is the local electron number density. Electrons that are not part of the relativistic fluid can be trapped in the wave, 
under certain conditions, and accelerated. 

It is hoped that electron beams from LWFAs will eventually provide sources of relativistic electron beams and x- 
rays for applications, such as chemical and biological spectroscopy, medical imaging and radiation therapy. If viable, 
they are expected to revolutionize such sources, making them available at lower cost to small scale facilities such as 
universities or hospitals. 

With advances in laser technologies, pulses of duration close to the plasma period (r p = 2tt/u p ) can now be generated 
directly. It was suggested that a high- amplitude wakefield could be generated by multiple suitably spaced short pulses 
[2]; if subsequent laser pulses are timed to push the electrons in their direction of motion as they overshoot their 
equilibrium position, then the plasma wave amplitude can be grown resonantly. However, as pulse length reduction 
has been accompanied by an increase in the power of the short pulses, multiple pulses have been deemed to be 
unnecessary since it is now possible to grow large amplitude plasma waves with a single intense ultrashort laser pulse. 

These high- amplitude waves can be a source for the trapped electrons by self-injection from the plasma itself. This 
occurs in three dimensions when the accelerating phase of the electric field reaches sufficient amplitude to prevent 
electrons from slipping into a decelerating phase of the wake [3] . In recent experiments [H |5] , it has been demonstrated 
that a single pulse injection process can produce electron beams of narrow energy spread and small divergence, but 
requires > 10 TW < 40 fs laser systems, which are by no means ubiquitous. In addition, since the propagation of a 
pulse — with IX 2 > 10 19 Wcm -2 required for self-injection — is dominated by modulational effects, a higher degree of 
complexity may be required to control the interaction. 

Alternative means of injecting electrons by the influence of secondary laser beams have also been previously sug- 
gested [H H IE] • These include a scheme whereby counter-propagating laser pulses have been shown to control some 
of the electron beam properties [9 . This demonstration has spurred further interest in optical injection in LWFAs. 

In this paper we present the first evidence for production of mononenergetic electron beams by all optical injection 
in a dual collinear ultrashort laser pulse geometry (i.e. with pulses overlapping, parallel and traveling in the same 
direction). We note that though collinear pulses were used in [lOj [TTj, in those the resulting spectrum was not 
monoenergetic and the mechanism was different. The resonant growth of the plasma wave is experimentally inferred 
by the production of monoenergetic electron beams by the two laser pulses, when correctly spaced, as neither of the 
pulses on their own do this. 

A laser pulse with a large focal spot (/-number) was used to generate a low amplitude plasma wave. This was then 
used as a guiding structure for a tightly focused (low /-number) pulse. Particle-in-cell simulations of the experiment 
show that the tightly focused pulse was guided at a smaller spot-size than the guiding pulse, for much longer than its 
Rayleigh range if correctly phased within the wakefield [12] . Electrons were produced when the low /-number pulse 
focuses, and are accelerated in the resonantly driven wake, resulting in monoenergetic beams. 

The experiments were carried out on the 600 mJ arm of the tl = 40±5 fs full-width half-maximum (FWHM) 
Ti: Sapphire Astra laser operating at wavelength Ao = 800 nm. A thin (5 mm) beam-splitter was used to produce 
two E pu i se = 300 mJ pulses that were focused using //3 and //16 off-axis parabolic mirrors collinearly onto the edge 
of helium flow from a 2 mm diameter supersonic gas nozzle (Fig.[T]). On ionisation by the leading edge of the laser, a 
plasma of electron number density n e = 1 x 10 19 cm -3 was produced. This density was chosen so that the pulse with 
//16 focusing did not produce a measurable electron signal, but so that electrons were detected using //3 focusing. 
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FIG. 1: (color online) Schematic of the experimental setup. 
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FIG. 2: (a-c) Quasi- monoenergetic electron spectra with highest charge shown for various timings between the drive and guiding 
laser pulses, compared with d) typical electron spectra with the drive pulse only. n e — 1 x 10 19 cm -3 and E pu i se < 300 mJ. 



The full width half maximum (FWHM) spot sizes for the pulses focused by the //3 and //16 optics were 5 jam and 
25 /im respectively. This resulted in focused intensities of 1.9 x 10 19 Wcm -2 and 0.8 x 10 18 Wcm -2 corresponding 
to normalized vector potentials ao of 3 and 0.6. In both cases, the longer focal length pulse passed through the 
beamsplitter, for which a 5 fs increase in pulse length was measured. A motorized timing slide allowed control of the 
longitudinal spacing of the two pulses. The pulses could be overlapped to an accuracy of about one focal spot size 
(i.e. ~ 25 /im), and were limited by the mechanical stability of the laser pointing at the time. 

The timing slide position was measured to ±1 /im (±3.3 fs) accuracy. Accuracy was limited by the measured 
relative timing of the two pulses At, which was done by maximizing the plasma defocusing of one pulse by the other, 
at a fraction of a percent of full power. The error in this method was therefore of the order 40 fs (i.e. ~ tl). The 
energy spectrum of the accelerated electrons was obtained using a magnetic spectrometer with image plates as the 
detector, as in [4]. The acceptance cone of the electron spectrometer was //200. Light emitted from the interaction 
was reimaged orthogonally to the beam axis, to discern the propagation of the laser. 

Using the pulse with //16 focusing — henceforth referred to as the guiding pulse — alone resulted in no measurable 
electron signal above noise level at the chosen density. The pulse with f/3 focusing — henceforth referred to as the drive 
pulse — would consistently produce electron bunches of a few pC charge, but these were consistently best described 
by a single temperature (i.e. one dimensional maxwellian momentum distribution) fit. 

Using both pulses, the spectra obtained were non-maxwellian with obvious structure (Fig. [2^i-c) . Indeed on many 
of the shots the maxwellian spectra that the drive pulse alone produced were almost completely suppressed, with 
the majority of recorded electrons in a single low energy spread beam. For example, the beam shown in Fig. [2J3 had 
an energy spread AE — 0.6 MeV (FWHM), which is the narrowest energy spread beam that we have yet measured 
from a self-injecting plasma accelerator. This monoenergetic beam was also very well collimated, with an angular 
divergence less than the opening angle of the collimator, 2 x 10 -5 steradians solid angle or 2.5 mrad half-opening 
angle. This also indicates that the measured energy spread in this case is also determined by the actual size of the 
beam, and so was probably less than quoted. 

The relative mean energies, Wr, and total charge in the spectrum are plotted for all shots in the data set in Fig. [3J 
including non- monoenergetic spectra. The mean energy W = J °° W dN(W)/ J °° dN(W). Wr = W — min(VK^). 
The total charge has a maximum for At = 30 fs, which is likely to be the timing when the pulses were overlapped 
sufficiently to act as a single pulse with double the energy, since it is within the error by which At = was set 
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FIG. 3: (color online) a) Relative mean energy Wr and b) total charge Qt as a function of the relative arrival time between 
the pulses for all shots. 



(Fig#). 

From Fig. [3^i, it is evident that there is also an enhancement in the mean energy, but for a larger range of timings. 
These represent timings from when the pulses are overlapped, to when the guiding pulse leads the driver by a few 
plasma periods (« 6). This indicates the extent of the coherent plasma wave structure for this density. Note that due 
to the low intensity of the guiding pulse, a single wave-period is not expected as in [13]. However, the decay of the 
energy enhancement indicates the eventual damping of the coherent plasma wave due to, e.g. radial anharmonicity 

The oscillatory behavior of the weighted energy with increased timing appears to have a well defined frequency. 
However this is likely to be an artifact. Periodic behavior is expected due to the underlying plasma wave structure, 
but the frequency here is too low at oo « uo p /2. However, the pulse length (~ 12 jim) was longer than X p ~ 7.5 /xm. 
Hence, it is likely that the periodic behavior is smeared out by the limiting pulse length. In addition, in our set-up 
the error in the relative timing was on the order of X p due to variations in pulse length and positioning variation. 
This is reflected in the scatter of data points in Fig. [3] The oscillatory behavior is probably a combination of these 
errors and also the aliasing of an underlying frequency (i.e. of the plasma wave) which is higher than that in the 
measurements. Nevertheless, it is worth emphasising that Wr can be enhanced but only when the driver follows 
within a time shorter than the decay time of the guide pulse wake. Yet over the same time-window, small variations 
in At can also lead to a decrease in Wr. 

Another difference for the interaction with two beams is the images of side-scattered light. It has been demonstrated 
that this side-scatter radiation is primarily Raman side-scatter, and its length is an indication of the initial evolution 
phase of the high intensity laser plasma interaction [15] . With the drive pulse alone, the extent of the emission is 
~ 10 /im, which is comparable to the Rayleigh range of this pulse. Side-scatter of the guiding pulse by itself was not 
measurable above noise level. However, when the drive pulse overlaps the guiding pulse, an extended emission region 
on the order of 1 mm is observed, which is longer than the Rayleigh range of the guiding pulse. 

An extensive series of two dimensional particle-in-cell simulations were run using the code OSIRIS [16]. These were 
run under similar parameters to the experiments and for a large number of relative arrival times, At a between the 
pulses. Due to computational constraints, At a could not be longer than a few plasma periods, however the control of 
the timing within that range was obviously much higher than in the experiment. 

The simulations clearly indicate the presence of plasma wave guiding [12] of the tightly focused pulse. This coun- 
teracts the filamentary behavior of the tightly focused pulse and extends its propagation as a single filament. There 
is a strong phase dependence to this effect in these simulations arising as a result of the fine control of At a . Plasma 
wave guiding relies on the guided pulse being in a predominantly focusing refractive index structure. Since a plasma 
wave is periodically focusing and defocusing it is evident that for certain At a the pulse will be guided and for a tt/2 
phase shift the pulse will be defocused. 

For simulations when the drive pulse is in phase with the plasma wave created by the guiding pulse, it is focused 
by the density depression and its propagation is dictated by the guiding pulse. In addition, the two wakefields are 
resonant, and therefore the combination a larger amplitude wake than either pulse alone. When it is out of phase 
with the wavebucket, the density structure acts to defocus the pulse. This guiding and defocusing is shown in Fig [4] 
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FIG. 4: (color online) The normalised intensity e 2 Eq / m 2 c 2 ujq of driver pulses having propagated 1500/ira in a background 
density of n e — 1 x 10 19 with a separation of a) 21/xra (well phased) and b) 18/xra (not well phased). Contours are given at 
intervals of 0.25. 
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FIG. 5: (color online) Electric field of pulse averaged longitudinally (z) and then plotted as a function of propagated distance 
for a) driver with no guide pulse, b) driver with guide pulse and c) the guide pulse. The graphs show (left) the intensity profiles 
in the transverse (x) direction and (right) the wakefield and pulse profiles along the axis of propagation, d) is the laser envelope 
and e) is 5n e /n e o for the wake. 



In Fig [4]3 the wake of the guiding pulse is out of phase with the drive pulse. The drive pulse is mainly defocused with 
only a small proportion of laser energy trapped in the density well. However, Fig [4^i shows the situation when the 
drive pulse phase is well matched within the guiding wake, and hence a large fraction of the pulse energy is trapped 
in a single stable filament. 

In Fig. [5j time histories of the intensity profiles of (a) the unguided drive pulse, (b) the guided drive pulse and (c) 
the guide pulse are shown as they propagate through the plasma. The guided drive pulse (b) is in phase with the 
wake, and is clearly guided over a distance longer than its Rayleigh length (zr = 33 jim). The non-ideal pulse lengths 
(similar to those of the experiment) meant that not all of the driver is within the focusing phase of the plasma wave, 
and hence it diffracts away. However, the majority of the pulse energy is trapped in a single filament. 
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The inset graphs to Fig. [5] show the intensity profiles for the pulses and wake density profile after 1.05 mm prop- 
agation. A gaussian fit to the temporally averaged central filaments gives a full width at 1/e 2 , 2w for each of the 
pulses. The guide pulse (c), is very close to X p in spotsize, 2w — 9.0 /im, as is expected for a self-guided short pulse 
[l5j[T7]- For the unguided driver pulse (a), the spotsize is 2w = 20.8 /im and has significant wings, whereas for the 
guided driver pulse (b), the spotsize is 2w = 5.6 fim. 

Earlier in the interaction (a and b, inset boxes), the drive pulse focuses to high intensity and traps an electron 
bunch. After defocusing, electron trapping ceases, and the behavior is strongly influenced by the presence of the guide 
pulse. In the unguided case, the lack of quasistatic wakefields maintained for a long propagation length results in a 
broadening of the electron spectra c.f. [15]. However, in the guided case, the electron bunch is accelerated in a uniform 
wakefield, resonantly generated by the two pulses for a relatively long time, and therefore remains monoenergetic. 
This means the electron trapping can be controlled through the focusing of the drive pulse. 

In summary, the production of monoenergetic electrons has been demonstrated for the first time by the interaction 
of two collinear pulses, through plasma wave guiding and resonant plasma wave growth. Under the right conditions 
of laser beam overlap and timing, electron beams of high quality, both in terms of energy spread and divergence, were 
produced. By controlling the propagation of a tightly focused driver pulse, this scheme offers the ability to control 
monoenergetic electron beam production and is easy to implement. Consequently this scheme may be of interest in 
the development of compact sources of energetic electron beams. 
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